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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE I 9 O 3 

EXPERIMENTAL INVESTIGATION OF MOVING PRESSURE 
DISTURBANCES AND SHOCK WAVES AND CORRELATION 
WITH ONE-DIMENSIONAL UNSTEADY-FLOW THEORY 
By Paul W. Huber, Cliff E. Fitton, Jr., and F. Delpino 


SUMMARY 


Experimental pressure— time measurements, along with schlieren 
photographs of the unsteady— flow phenomena initiated in constant— cross- 
sectional— area passages, are presented. The methods of calculation of 
the flow phenomena are based on the assumptions of one-dimensional, 
nonviscous, unsteady flow. The unsteady flow — which is studied both 
experimentally and theoretically — includes propagation of pressure 
disturbances, which may be either expansion zones or discontinous 
compression fronts, temperature contact discontinuities in the flow, 
and the reflections and interactions of these flow quantities. A very 
close agreement is shown between the experimental measurements of the 
flow and the flow predicted by the theory over a wide range of flow 
conditions. Special instrumentation was developed during these tests 
which greatly facilitated the measurement of the unsteady-flow phenomena. 


INTRODUCTION 


The important role played by moving pressure disturbances in 
various aerodynamic conf igurations is becoming increasingly apparent. 

In pulse— jet engines and heat engines, such as the comprex or pressure 
exchanger (reference l), moving compression and expansion waves are 
alternately used to compress the gas before heat is added and to expand 
the hot gases in the same passage, in much the same manner as is obtained 
by the action of a mechanical piston. The action of pressure disturbances 
in the subsonic part of a supersonic diffuser is recognized as causing an 
instability of the shock wave with a resulting adverse effect on pressure- 
recovery efficiency. In some cases (reference 2), constant— geometry 
supersonic diffusers have been started at Mach numbers below the theoreti- 
cal starting Mach number (based on steady— flow considerations) by means 
of pressure disturbances. In steady— flow— type burners, as well as in 
combustion engines, the action of pressure disturbances and detonation 
waves produce a marked change in the combustion phenomena. 
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Although these actions of pressure disturbances in aerodynamic 
and thermodynamic phenomena have become recognized and their effects 
have been partly evaluated from the general trend of experimental tests, 
the amount of available experimental data has been insufficient to 
verify the theory of unsteady flow. Tills lack of fundamental experi- 
mental data has been partly due to a lack of instrumentation capable of 
responding with any degree of accuracy to the rapidly changing conditions 
of this type of flow. 

In general, many of the one— dimensional unsteady— flow cases of 
fundamental interest can be treated mathematically with existing methods. 
(See references 3 to 6.) These methods, in their present form, however, 
are not applicable to viscous flows, flows involving heat transmission, 
or flows of more than one dimension, which are usually the flows of 
engineering interest. 

The purpose of the present investigation is to compare the results 
of an experimental study of one— dimensional unsteady flow with results 
based on the theory of one-dimensional, nonviscous, unsteady flow. A 
better knowledge of the gas dynamics of this simple type of unsteady 
flow should prove to be of great assistance in the study of the more 
complicated flows. As a secondary objective, correlation of theory and 
experiment has been used to determine the value of instrumentation for 
measurements of one— dimensional, nonviscous, unsteady flow. 


SYMBOLS 


P 


pressure, atmospheres 


u 


flow velocity, feet per second 


a 


velocity of sound, feet per second 


t 


time, seconds 


T 


7 


temperature, degrees Fahrenheit 

ratio of specific heats at constant pressure to specific heats 


at constant volume 


x 


distance, inches 


M 


Mach number in stationary shock 



Y 


flow velocity in stationary shock coordinates, feet per second 


W 


velocity of propagation of shock, feet per second 
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E voltage, volts 

l length., inches 

Ap pressure change at crystal due to given pressure front 

Subscripts: 

a standard atmospheric (29.92 in. Hg) 

0 high— pressure chamber before diaphragm burst 

1 low-pressure chamber before diaphragm burst 

b before reflection 

c after reflection 

22 between temperature contact discontinuity and expansion, 

before reflection 

2 between temperature contact discontinuity and shock, before 

reflection 

3 between reflected shock and end of low-pressure chamber 

4 between reflected expansion and end of high-pressure chamber 

L low-pressure chamber 

H high-pressure chamber 

Superscripts: 

1 high-pressure side of shock 

" low-pressure side of shock 


THEORETICAL METHOD 


The fundamental relationships for the propagation of plane 
disturbances of finite amplitude were first derived by Riemann in 1859 
for one— dimensional isentropic flow. (See references 3 an-cL 7 .) The 
derivation showed that the propagation of a large disturbance can be 
treated as the propagation of an infinite number of small disturbances 
where each small disturbance is propagated at a velocity with respect 
to the fluid equal to the velocity of sound at that point. If the 
fluid at that point has a flow velocity of u, then the velocity of 
propagation of this point on the disturbance is a ±. u, where a is 
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the Telocity of sound and the plus sign or minus sign is used if the 
point is being propagated in the direction of flow or against the 
direction of flow, respectively. For a streamline of constant cross- 
sectional area and for a fluid which has a constant value of 7 , a 
disturbance propagated in the direction of flow was shown to produce 


no change in the quantity of the flow 



through which it 


has been propagated; whereas a disturbance propagated against the direc— 

2 

tion of flow was shown to produce no change in the quantity u + — — a. 

The propagation of expansion fronts on these disturbances was shown to 
be such as to make the expansion fronts less steep with time (that is, 
spreading expansion zones); whereas compression fronts were shown to 
become more steep with time and eventually to become discontinuous 
pressure fronts (shocks). From this result, it can be concluded that 
the shape (the shape being the pressure distribution along the axis of 
propagation) of any finite isentropic disturbance continually changes 
with time. In order to calculate the shape of a disturbance at any 
given point in time, therefore, the shape must obviously be known at 
some earlier time. If an experimental measurement on the moving 
disturbance is not trusted, then the time and conditions under which 
the disturbance was generated must be known (time = 0) . 


Classical methods for the treatment of disturbances from initial 
conditions which are known include the following: 

1. Assumption that a piston rapidly accelerates to a given velocity 
which is then held constant. This acceleration of the piston generates 
a compression front which moves ahead of the piston, while the fluid 
behind the front flows at the velocity of the piston. Deceleration of 
the piston to rest consequently generates an expansion zone which moves 
ahead of the piston and reduces the fluid velocity to zero. 

2. Assumption that a diaphragm which separates a chamber of high 
pressure from a chamber of low pressure is suddenly burst. In this 
treatment, a compression front is propagated into the low-pressure 
chamber; whereas an expansion zone is propagated into the high-pressure 
chamber. The movement of the compression front is such as to compress 
the fluid in the low-pressure chamber to an intermediate pressure and 
to accelerate the compressed fluid in the direction of the front. The 
movement of the expansion zone is such as to reduce the pressure in the 
high-pressure chamber to the same intermediate pressure and to accelerate 
the expanded fluid in a direction opposite to the front. A temperature 
contact discontinuity appears behind the compression and moves in the 
direction of fluid flow. This surface of discontinuity is the boundary 
between gas particles expanded from the high-pressure chamber and particles 
compressed from the low-pressure chamber. 

The diaphragm method was chosen because it seemed somewhat easier to 
reproduce experimentally. 
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In 1932, Schardin (reference 8) performed calculations from which 
the pressure distribution along tubes of constant cross section can be 
found at certain times after the shattering of the diaphragm. These 
calculations are based on the assumptions that the flow is nonviscous 
and that a plane shock wave is immediately formed which moves into the 
tube of low pressure while an expansion wave moves into the tube of high 
pressure. Schardin used the Rankine-Hugoniot relations for the shock 
fronts and Riemann’s relations for the expansion fronts. 

By an extension of Schardin’ s calculations to include reflections 
and interactions of the waves, as well as the effect of the temperature 
discontinuity which exists in the flow, a complete time history of the 
phenomena can be calculated. In order to find the conditions which exist 
immediately alter the shattering of the diaphragm, the method used by 
Schardin should be briefly reviewed. 

In figure 1 is shown the pressure distribution at typical times 
along constant-area passages which were separated by the diaphragm at 
t = 0. Figure l(a) shows the pressure distribution along a tube in 
which a diaphragm separates a pressure p^ from a lower pressure p-^. 

Figure l(b) shows the pressure distribution along the tube at a time t-j-, 

after the diaphragm has been shattered. A shock front is assumed to be 
propagated towards the region of low pressure and to be compressing the 
gas at pj to an intermediate pressure pg. At the same time, a zone 

of expansion is assumed to be propagated towards the direction of high 
pressure and to be reducing the pressure Pq to the same intermediate 

pressure pg. In order to find all the flow quantities existing in the 
tube at the time t-^, the following assumptions are made: 

(1) The value of 7 is constant. 

(2) Velocities are positive in the positive x-direction. 

(3) The tube is of constant cross-sectional area. 

(4) Before the diaphragm is shattered, the temperature in the 

tube is everywhere the same and the velocity is everywhere 
equal to 0. 

(5) The flow is nonviscous and involves no transmission of heat. 

Immediately after the rupturing of the diaphragm, a temperature 
contact discontinuity exists in the flow somewhere between the compression 
front and the expansion front as mentioned previously. This discontinuity 
also represents an entropy discontinuity, because this discontinuity 
represents no pressure increment and, hence, no velocity increment. In 
order to satisfy the assumption that an expansion zone is moving in the 
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negative x-direction, it now follows from Biemann’s work that thfe 
velocity U2 must he in the positive x-direction. Therefore, 


u 0 + 



Ug + 



a 22 


which may he expressed as 



( 1 ) 


because the process is isentropic and 

7-1 

a 22 _ /PA 27 

a 0 V P 0/ 

Equation (l) expresses the velocity increment across the expansion zone 
as a function of the pressure ratio across the zone and aQ, which is a 
function of the initial temperature. 

The velocity increment across the shock front can he expressed in 
terms of the pressure ratio across the front and a b by using the 

Bankine— Hugoniot relations for a stationary normal shock as follows. 

The coordinate system in figure l(b) can he changed to reduce the 
shock velocity W b = 0 by expressing the gas velocities on either side 

of the shock as 


V" 


= U 1 " W b 


and 


V* 


u 2 - W b 
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where V" is the gas velocity on the low-pressure side of a stationary 
normal shock and V* is that on the high-pressure side. (See fig. 1(c).) 
The velocity increment across the shock is then 


V» - V" = -M 2 a 2 “ ( -M l a l') 


which may he written as 



(2) 


where and Mg are the Mach numbers on the low-pressure and high- 
pressure side, respectively, of a stationary normal shock and from the 
Kankine-Hugoniot relations. 



Also, from the stationary— shock relations. 


M n = 


(7 + l) — + (7-1) 
?! 

27 


and 


Mg = 


(7 + 1)^ + (7 - 1) 

? 2 

27 
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so that equation ( 2 ) may he -written a 3 



Equation (3) is the expression for the velocity increment across the 
shock front, for either coordinate system, as a function of the pressure 
ratio across the front and a^* It might he pointed out that the selection 
of coordinate systems other than the one used herein yields the same 
relation. (For example, see reference 5.) Equation ( 3 ) may he equated 
to equation (l) because a^ = aQ and u^ = uq = 0 as follows: 


Pi 

Po 



— - 1 
Pi 


- f7 

7-1 


27 


(y + 1 22. + 1 


'7 - 1\7 - 1 Pq 


(4) 


Equation (4) can easily he solved graphically for P 2 /P 3 / (See 
fig. 2.) By substituting pg/p^ in equation ( 3 ), Ug may he found 
because is known from the initial temperature in the tube. Also, 

may he found from 





or the pressure ratio across the shock may be written as a function of 
the shock speed 



7-1 


when uj = 0. 
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All the conditions in the tube are now known for any time tp after 
the shattering of the diaphragm until reflections occur at the ends of 
the tube. These conditions satisfy the original assumptions of an 
expansion zone moving in the negative x— direction and a compression shock 
moving in the positive x— direction, as well as a temperature contact 
discontinuity (&22 4 & 2 ) moving at the velocity U2 across which there 
is no change in pressure or flow velocity. The pressure distribution 
along the tube is easily found because the shock has moved a distance Wptp 
and the expansion zone AB in figure l(b) has moved distances of -&o^b 
and — (a 2 2 — UgHp for points A and B, respectively. The expansion 

zone becomes less steep with time because point A moves faster than 
point B; whereas a compression shock exists immediately and moves in the 
positive x-direction because the propagation velocity behind the 
shock (up + a^) is greater than that ahead (aq) „ 

Figure l(d) shows the pressure distribution along the tube at a 
time t c after the shock has been reflected from the low-pressure end of 
the tube and the expansion has been reflected from the other end. The 
conditions existing after reflection of the shock can be found in a 
manner similar to that used to obtain equation (3); that is, the reflected 
shock velocity — W c = 0 can be reduced so that the flow velocities on 
either side of the shock become (fig. 1 (e)) 


The velocity increment across the shock can therefore be expressed in 
terms of its pressure ratio and ag as follows: 


V" = u 2 - (-W c ) 


and 


v» = u 3 - (-W 0 ) 





( 6 ) 
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Because U 3 = 0 at the end of the tube, P 3 /P 2 can expressed 
as a function of pg/p-^ by combining equation (6) with equation ( 3 ) where 


&2 

a l 



The resulting expression is, for 


7 = lA, 


£3 

p 2 


8^-1 
_Pl 

6 

P 1 


( 7 ) 


The shock velocity W c can be found from 

— W c = u^ — Y* = M^a^ 
or 

-¥ c = "2 “ Y " = 


where and are the Mach numbers on the high-pressure and low- 

pressure side, respectively, of the reflected shock in stationary 
coordinates; that is. 


t 


M- 


l z + 1 £2 + 1 
7 - 1 p 3 
27 

7-1 


= 


7 + 1 £3 

[7 - 1 pp 

27 


+ 1 


7-1 
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and 




+ 1 
7-1 


P3 

P 2 


7+1 

7-1 


All the flow conditions after the reflection of the shock are now 
known. Schardin’s calculations (reference 8) are extended to the reflec- 
tion of the shock. (See table I.) 

In order to extend the calculations to a complete history of the 
flow inside the tube, the next step is to find the reflection of the 
expansion zone from the end of the tube. The conditions after total 
reflection of the expansion are easily found from Riemann’s equation 
where 


2 2 
- — ~ a b = u 2 “ — X a 22 

which may be solved for a[^ because u^ = 0 at the end of the tube. 
The pressure p^ can then be found from the isentropic relationship 

27 
7^ L 


Because of its thickness, however, the total reflection of the 
expansion requires a finite time so that the pressure history at the 
end of the tube during reflection, or the pressure distribution along 
the axis of propagation between points A and B must be found somewhat 
differently. This pressure history is obtained by considering the 
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expansion zone AB to consist of an arbitrary number of small discontinuous 
isentropic expansion fronts (the usual characteristics method) rather 
than to be a smooth pressure variation between zone AB. For an illustra- 
tion of this method, see figure 3 in which the expansion zone AB is 
assumed to consist of two expansion zones Y and Z which are replaced with 
isentropic fronts Y and Z. This method is, of course, not rigorous because 
Kiemann has shown that an expansion front cannot remain discontinuous; 
however, by the choice of small steps the method provides a good approxima- 
tion. Each front can now be considered individually, and the conditions 
of u and a on either side of each front can be found from the initial 
pressure increment assigned the front by 


u 


+ 2 
"7-1 


a 


= Constant 


* 


In figure 3, step 2, p = 1.70 atmospheres was arbitrarily assigned 
the region between fronts Y and Z. From the relation u + 5a = Constant 
across either front Y or Z (because the fronts are being propagated 
against the fluid flow) and the adiabatic relations, the value of u 
and a in this region are found. The speed of each front is assumed to 
be an average of the values of u ±. a on either side of the front. For 
this reason, the characteristics method represents an approximation. 

The total reflection of the expansion from the end of the tube is 
obtained by reflection of the small discontinuous expansion fronts 
individually and by consideration of the subsequent interactions of the 
reflected with unreflected fronts. In figure 3* step 3, the conditions 
at the tube end after reflection of front Y are first found by using the 
relation u — 5a across reflected front Y and by knowing that u = 0 at 
the tube end. The conditions between fronts Y and Z after interaction 
are next found by using the relationships u + 5a = Constant across 
front Z and u — 5a = Constant across front Y and by obtaining a 
simultaneous solution for u and a between fronts after the interaction. 
(See fig. 3, step *<-.) The final step is the reflection of front Z 
from the tube end. (See fig. 3, step 5*) 

Because the speeds of all the fronts are known and because all the 
fronts originate at the diaphragm at t = 0 , a more convenient method of 
illustrating the positions of the various fronts in the tube at any time 
is to plot the time histories as a diagram of x against t as in 
figure k where the expansion is divided into seven fronts. Pressure 
variations, as a function of time or distance along the tube, are easily 
taken from such a plot because conditions in each enclosed area are 
known and are constant throughout the area. In figure 5j the preceding 
method was used to show the pressure distribution along the tube between 
points A and B at a time t^,. In figure 6 , this method was used to show 
the variation of pressure at the end of the tube with respect to time at 
reflection of expansion zone. It can be seen from figures 5 and 6 that 
the pressure variation inside the expansion front is nonlinear. 
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Figure 4 shows that, depending on the ratio of the lengths 

the temperature contact discontinuity is overtaken "by either the reflected 
expansion zone or the reflected shock front.' Whichever is the case, it 
can he shown that the pressure ratio of the front is altered after passing 
through the temperature contact discontinuity which, in turn, generates 
a second front that propagates itself in the opposite direction. The 
method used to obtain the conditions which exist after the front has 
undergone this interaction with the temperature contact discontinuity is 
similar to the method already used. First, the assumption is made that 
a second front is generated and travels in the opposite direction from 
the original front. Then, knowing that there is no velocity increment 
or pressure increment across the new discontinuity necessitates the 
finding of four unknown conditions: namely, the pressure and velocity 

existing in the region between the two fronts, and the two temperatures. 

These conditions can be found in the same manner in which the condi- 
tions after the shattering of the diaphragm were found, except that the 
generated front must be assumed to be either a shock or an expansion. 

After solutions are obtained and the assumption is indicated to be 
incorrect, a second solution must be obtained,, If the interaction under 
consideration is that of an expansion zone with a discontinuity, the 
expansion zone must be considered as a number of small fronts (as was 
the treatment for the reflection of an expansion zone from the end of 
the tube) in order to obtain a pressure variation with time or distance 
inside the expansion zone during the interaction. This method is rather 
laborious and is not required when solving for conditions after the 
interaction but should be used for an accurate history of events in a 
diagram of x against t. Interactions of expansion zones and shock 
fronts with discontinuities have been analyzed by this method, and, in 
general, the generated front is very weak in proportion to the original 
front. An interaction of a shock front with a discontinuity is shown in 
figure 4 in which the generated front is a compression traveling in the 
opposite direction from the original shock. 

In order to be able to find conditions in the tube at any time or 
position, the one remaining phenomenon which must be considered is an 
interaction of a shock front with an expansion zone. For this case, 
the conditions which exist after the interaction is over are again 
readily found by the methods used for finding the conditions after the 
shattering of the diaphragm. The only assumption that satisfies a 
solution of this interaction is one in which a shock and an expansion 
are propagated in the same direction after the interaction as before 
the interaction and in which the flow between the two fronts is reversed 
in direction. The values of u, a, and p, which exist between the 
fronts after interaction, can be obtained by making this assumption and 
by obtaining simultaneous solutions of the shock equations and the 
i s entr op ic— front equations. The expansion, however, must be considered 
as a number of small fronts and the method must be applied the same 
number of times if a pressure variation of time or distance is to be 
obtained during the interaction. An interaction between a shock front 
and an expansion zone is shown in figure 4 . 
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By use of the methods which have heen discussed it is possible to 
determine al 1 the subsequent phenomena and to obtain the complete pressure 
history in the tube at any point in the tube or a complete pressure 
distribution along the tube at any time after the diaphragm is shattered. 
The assumptions of nonviscous flow of a gas which has a constant value 
of 7 and which involves no heat transmission should be valid in many 
cases. However, for passages of very small cross section or long length 
and for pressure variations involving large temperature variations, 
these assumptions would undoubtedly lead to incorrect results. The 
method which has been presented can, in many cases, be greatly simplified 
for small pressure variations. In these cases the shock can be treated 
as an isentropic pressure front, and Riemann's relations, rather than the 
shock relations, can be used with very little error. Also, the expansion 
zone need not be broken into a large number of weaker fronts; usually two 
or three will suffice. 


EXPERIMENTAL METHOD 


A diagrammatic sketch and photographs of the apparatus are shown in 
figures 7 and 8, respectively. A description of the design and function 
of the instrumentation, much of which was developed specifically for 
these tests by the Langley Instrument Research Division, was prepared 
by Mr. F. Delpino and is given in the appendix. 

The apparatus in which the waves were produced consisted of four 
passages of constant internal cross-sectional area which could be 
connected together in any sequence by means of clamps and by mating 
flanges to form one long passage with smooth concentric inner surfaces. 
The internal dimensions were everywhere 3 inches wide by 3 inches high. 

3 

The walls were fabricated from jg--inch brass stock by brazing the joints 

along the outside surface of the passage. The four sections included a 
section 8 inches long which contained a solenoid— actuated trigger used 
to puncture the diaphragm. (See fig. 7(c).) Also included was a 2k — inch 

section, a 48-inch section, and a 24^-- inch section, which was equipped 

with two pairs of glass windows of such dimensions that the full 3 —: inch 
height of the passage for a distance of 6 inches along the passage 
could be photographed from either pair of windows. The glass windows 
were mounted in such a maimer that the inner surfaces of the glass were 
flush with the inner surfaces of the passage with gaskets between the 
joints to prevent leakage of air. Each section of passage was equipped 
with two amal 1 static-pressure orifices (0.04— inch diameter), one which 
could be connected to a pressure gage and the other, to a source of 
pressure or vacuum. At the positions shown in figure 7(a), each section 
was equipped with tapped holes along the top of the passage which could 
accomodate either a blank plug screwed in flush with the inner surface 
or a plug in which the pressure-measuring instrument was mounted flush 
with the surface as shown in figure 7(b). Three metal end plates were 
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provided of thiclaiesses of 3/16 inch, 3/8 inch, and ll/l 6 inch, each of 
which could he clamped to the end of any section of passage. 


In all the tests, the 8 — inch trigger section, the 24— inch section, 

3 

and the 5 - — inch end plate were clamped together to form what is called 
o 

the high-pressure chamber e The low-pressure chamber then consisted of 

either or both of the remaining two sections in addition to the 11 — inch 

lo 

end plate. In this manner three different lengths of low-pressure chamber 
were available in four possible combinations of section — that is, one 


24-1- inch, one 4-8 inch, and two 72 -- inch 


in which the glass section* 


could be at either end. All the clamped joints between sections and 
between a section and an end plate were equipped with rubber gaskets of 
such dimensions that the inner edges were flush with the inner surfaces 
of the passages. The trigger was equipped with a sliding pressure seal 
and a soft rubber cushion to damp out mechanical vibration. (See fig. ?•) 
The blank plugs and instrument plugs were equipped with rubber gaskets. 

The diaphragm used to separate the high-pressure chamber from the low- 
pressure chamber consisted of from one to four thiclaiesses of a special 
paper, depending on the pressure difference between the two chambers, 
and was inserted in this joint along with a rubber gasket before clamping. 
The particular paper used, American Blueprint Co. Inc., No. 4-5D, Type 2, 
tracing paper, was selected after many preliminary trials because of the 
following favorable characteristics: very little stretching under 

tension; minimum porosity and hence minimum leakage of air through the 
paper before puncturing; and, finally, disintegration into small pieces 
when punctured rather than simply tearing into four pie— shaped pieces. 

The pressure-measuring instrument was a commercial unit consisting of 
an ammonium dihydrogen phosphate piezoelectric crystal mounted in a metal 

case of — inch diameter and 1— inch length. The electrical characteristics 

8 

of this instrument are described in the appendix. This instrument has 
the following desirable characteristics for this type of test: rugged 

construction, high acoustic impedance, high sensitivity, and no cavity 
resonance. 


Before each series of tests, the high-pressure and low-pressure 
chambers were clamped together without a diaphragm and the entire unit 
was checked for air leakage. Before each test, the small bits of paper 
left in the chambers from the preceding test were blown out with an air 
hose. After a new diaphragm was clamped between the high-pressure and 
low-pressure chambers, the pressure in each chamber was adjusted to the 
desired values and read from calibrated pressure gages to an accuracy 
of 0.1 pound per square inch. The pressures were held steady for a few 
minutes before each test so that the temperatures could come into 
equilibrium with the room temperature. Because the gage leads and 
pressure— supply leads were connected to the apparatus through small 
orifices and because the duration of the test was only a few milliseconds. 
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the small amount of flow Induced through these orifices during the test 
were believed to be of no consequence. 

Photographs of the flow were taken with a conventional schlieren 
system which included two parabolic front— surfaced mirrors of 6— inch 
diameter and a high— voltage— spark light source which was synchronized 
with the diaphragm burst as described in the appendix. Measurements of 
the absolute magnitude of the shock speed before its reflection from 
the end of the low-pressure chamber were made as described in the appendix. 

Preliminary tests were made by using a setup of smaller dimensions 
than the one previously described. These tests were made principally in 
order to develop the instrumentation. Additional preliminary tests were 
made by using the apparatus described in order to determine the effects 
of the following objectionable conditions on the pressure records: 

(1) Air leakage or air cavity around the crystal unit 

( 2 ) Mounting of the crystal unit in rubber instead of the Micarta 
mount shown in figure 7(i>) 

( 3 ) Mechanical vibration caused by the solenoid-actuated trigger 

(4) Failure of the crystal unit to be perfectly flush with the 
inner surface of the passage or end plate 

( 5 ) Failure of the diaphragm to disintegrate after bursting 

(6) Mechanical vibration of the end plates set up by the reflection 
of the shock waves 

Before the effect of these conditions on the test results are 
discussed, the method of interpreting the test records should be discussed. 
As pointed out in the appendix, the test records are essentially a plot 
of the voltage developed by the crystal unit as a function of time 
provided that proper electrical instrumentation is used. If, therefore, 
the pressure sensitivity of the crystal unit is known (that is, the 
voltage output as a function of pressure change at the crystal), the 
test results can then be compared with theory by comparing the theoreti- 
cal and experimental, pressure— time plots. The theoretical pressure— time 
plot for any position of the crystal unit can be taken from plots such 
as those shown in figures 9(a) to 9(c), which were obtained for a pressure 
ratio Pq/pi of 2.0 for the three lengths of low-pressure chamber, by 

using the methods outlined previously and by assuming the fluid to be 
air (7 = 1.4). In figure 10, the theoretical pressure— time plots are 
shown as dash lines for a number of different locations and lengths of 
low-pressure chamber, all for Po/Pp of 2.0. 

With the theoretical pressure— time plot known, the effects of the 
objectionable conditions previously listed were experimentally studied 
and the results of these preliminary tests are as follows: 
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(1) Records were taken with and without air leakage and an air 
cavity around the crystal unit. Comparison of the records indicated that 
a sharp fall of the trace occurred Immediately following a pressure rise 
at the crystal when leakage and air cavity were present. 

(2) Records were taken with the crystal unit mounted in rubier and 
were compared with records taken with the unit mounted in the Micarta 
mount shown previously. Oscillations of the trace on the records taken 
with the rubier mount occurred following a sharp compression. These 
oscillations indicated that the crystal was oscillating mechanically at 
the frequency of the rubier mount. 

(3) Records were taken with no pressure difference across the diaphragm, 
and a comparison was made of the trace obtained with and without a soft 
rubber cushion on the trigger mechanism shown previously. The mechanical 
vibration caused by the trigger hitting the stop on the side of the 
passage was found to be substantially eliminated by the rubber cushion. 

(4) Records taken with the surface of the crystal unit not flush with 
the inner surface of the passage or end plate indicated that the trace 
was distinctly altered because of the waves reflected from the crystal 
edges in this condition. 

(5) Records were taken of bursts in which more thicknesses of paper 
diaphragm were used than necessary, in which case the burst was more a 
tearing of the paper than a disintegration of the paper into tiny pieces. 
These results indicated that in cases where the paper did not disintegrate 
the rise of the trace as a result of a compression was more gradual and 
was followed by oscillations of the trace, probably a result of small 
waves behind the compression front caused by poor air flow past the 
remaining part of the diaphragm. 

(6) On a3 1 records taken with the crystal unit in the end plates, a 
rather violent oscillation of the trace would always occur following a 
compression. This oscillation would not occur follbwing an expansion, 
nor would it occur for either case if the crystal unit was mounted in the 
side of the passage. Also, the use of thicker end plates increased the 
frequency of this oscillation but had little effect on its amplitude. 

The addition of a 6-pound weight onto the plug in which the crystal was 
mounted (fig. 7 (b)) would decrease the frequency of oscillation for a 
given end— plate thickness. In all cases the rise of the trace due to 
the compression was reproduced (for a given pressure) if measured to the 
mean of this oscillation. The oscillation was a result of mechanical 
vibration of the end plate excited by the sharp pressure change in the 
compression and was extremely difficult to eliminate. 

The qualitative nature of all these effects can be predicted from 
general unsteady— flow theory; however, it was desired to determine the 
order of magnitude of the effects before analyzing the records. As a 
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result of these preliminary tests, the standard configuration for all 
subsequent tests, unless otherwise noted, incorporated a Micarta crystal 
unit mount as shown in figure 7(h), a soft rubber cushion on the trigger 

as shown in figure 7(c), an ^ — inch— thick end plate on the low-pressure 

chamber, and — as pointed out in the appendix — a 100-mi cromi cr of ar ad 
condenser across the crystal as well as a ^7-kilocycle filter in the 
preamplifier circuit. 


RESULTS AND DISCUSSION 


A large number of test records were taken at various crystal positions 
and with the three different low— pressure— chamber lengths. Records were 
taken over a range of pressure ratio Pq/Pi across the diaphragm from 1.3 

to 2.5 where p^, the pressure in the low-pressure chamber, was equal 

to 1 atmosphere or l/2 atmosphere. Typical pressure records which were 
obtained are shown in figure 11. 

Before the experimental records can be superimposed on a given 
pressure— time scale for comparison with theory, it is necessary to know 
the voltage output of the crystal unit as a function of the pressure 
change experienced by the crystal unit. This pressure sensitivity of 
the crystal unit was found experimentally by plotting all the voltage 
changes as measured from the test records as a function of the pressure 
changes experienced by the crystal as predicted by the theories previously 
outlined. In figure 12 and table II are shown the results of these 
voltage measurements, along with the positions on the test records at 
which voltage measurements were taken. These data indicate that, with 
only a very small scatter of points, a straight line can be drawn through 
the points starting from zero; that is, the output of the crystal unit 
is a linear function of the pressure change at the crystal, within the 
range of pressures shown. As shown in figure 13, the voltage measure- 
ments E were taken from the test records by first drawing a straight 
line (shown as a dash line) through the part of the trace immediately 
following either an expansion or compression and extending it to the point 
in time at which the voltage change began. This line corrects the 
measurement for the slight fall of the trace towards zero voltage inherent 
in the crystal and electrical circuit as pointed out in the appendix. 

No significant difference was noted among records taken with a number of 
crystals all of the same model type. 

Because the abcissa in figure 12 is the theoretical pressure 
difference experienced by the crystal unit for a given Pq/Pi across 

the diaphragm, it was necessary to determine experimentally if the 
absolute magnitude of this pressure difference was that predicted by 
theory. In order to do this without making a pressure measurement, the 
value of the shock speed W b was experimentally measured (see appendix) 
and the pressure p 2 was calculated from equation (5), p 1 and a^ 
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being known. The results of these measurements are shown in figure 2, 
where it can be seen from the measured values of pg/p^ that the actual 

values of Ap are only a few percent lower than the theoretical, values. 
The use of Ap to calculate the shock speeds gave values only a fraction 
of a percent lower than theoretical shock speeds. 


The experimental pressure— time records for 


&-2.0 
pi 


were super- 


imposed on the theoretical pre s stir e— time plots and are shown as solid lines 
in figure 10. In figure 10 the pressure scales (that is, voltage scales) 
for the theoretical trace were taken from figure 12, and the test records 
were superimposed by tracing directly from the records on a light table. 

The time scales and voltage against trace deflection scales were taken 
from the timing frequency traces and the alternating— current voltage traces 
on the test records. 


Figure 10 shows that the agreement between the theoretical and 
experimental pressure— time plots is very close. The slight departure 
which is shown is seen to be mainly due to the fall of the trace towards 
zero voltage, previously mentioned as a condition inherent in the electrical 

circuit. This fall was measured (see appendix) and was about 1^- percent 
of the voltage per millisecond. Although this correction was not made in 
figure 10, the agreement would be still closer if the correction were made. 

Schlieren photographs taken by using the 24^- — inch low-pressure 
chamber and by using the 72^- — inch low-pressure chamber with the glass 

section at the end plate are shown in figure 14. As a result of the 
schlieren photographs, shock— speed measurements, and pressure records, a 
number of interesting comparisons can be made between the experimental 
results and the theoretical predictions which were based on one— dimensional 
unsteady flow. 

The first assumptions made were that the flow was frictionless and 
that a plane compression shock is immediately formed after the shattering 
of the diaphragm. In figure 12, points taken at crystal locations 

ranging from 75- inches from the diaphragm up to 72^- inches from the 

diaphragm indicate that the shock does not exhibit any measurable 
attenuation within the scatter of the points (less than 1 percent of Ap). 
This assumption is further confirmed by the pressure— time records in 
figure 10 and the shock— speed measurements in figure 2 where attenuation 
would show up as a change in shock speed and thereby cause disagreement 
with the theoretical pressure-time plot. The schlieren photographs and 
pressure— time records show that the compression is not immediately one 
dimensional and discontinuous but that it becomes so very rapidly. 

At 7g inches from the diaphragm the schlieren photographs and pressure- 

time records (figs. l4(a), l4(b), and 10(d)) indicate that the compres- 
sion has a definite thickness; however, the schlieren photographs and 
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pressure— time records taken at 17 inches and 2h|- inches, respectively, 

(figs. l4(d), 10(a), 10(h), and 10 (c)) indicate that the shock is very 
sharp after this length of travel. This rapid steepening is to he 
expected from theoretical considerations because for this Pq/pi the 

value of u + a behind the shock is appreciably larger than that ahead 
of the shock (1462 fps as compared with 1126 fps). It cannot be definitely 
established from these results why the shock is not immediately plane 
and discontinuous, but it might be due to either an initial shock curva- 
ture because of the bowed diaphragm (in which case there would be oblique 
reflections from the walls trailing the shock) or due to the finite 
time required for the shattered pieces of diaphragm to spread out from 
the point of rupture. 

The next assumption made was that the temperature contact dis- 
continuity represented no pressure (that is, velocity) increment. The 
record in figure 10(d) shows that no pressure change is recorded at the 
point (about 1.75 milliseconds after the compression) where the dis- 
continuity passes the crystal. 

Finally, the theoretical method is based on the assumption that the 
diaphragm instantaneously disappears such that not only are plane discon- 
tinuous compression and expansion fronts generated but a plane temperature 
contact discontinuity is also generated. The schlieren photographs and 
pressure— time records, however, show that this assumption is not the case 
experimentally. The pressure— time records show that the expansion zones 
are slightly thicker than expected from theory (fig. 10 ), and the 
schlieren photographs and pressure— time records show that the temperature 
contact discontinuities appear to be a thick region of mixing in the flow. 
Figure lh, in which the reflected shock is being propagated through the 
temperature contact discontinuity, shows that some distortion of the 
shock is evident, probably a result of the unequal temperature distribution. 
From figure 10(a), a gradual rise occurs where a sharp compression was 
expected theoretically as a result of the interaction of this reflected 
shock with the temperature contact discontinuity. 

A possible explanation for these results might be that in the case 
of an expansion which is not initially plane because of the diaphragm 
bowing, oblique reflections of the expansion from the wall would, not 
theoretically catch up with the expansion. Also the temperature contact 
discontinuity, when not initially plane, would become a vortex in the 
flow because of its entropy variation. 

The pressure records have previously shown that the reflection of 
the shock from the end plates gave rise to mechanical oscillation of the 
end plates. The pressure waves generated as a result of this oscillation 
can be seen in figure 14 to be trailing the reflected shocks. 
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CONCLUDING EEMAKKS 


Experimental pressure— time measurements, shock— speed measurements, 
and schlieren photographs have "been presented and the results indicate 
that the unsteady flow generated experimentally conforms accurately with 
the flow defined by one-dimensional, nonviscous, unsteady— flow theory. 

The agreement of experimental measurements with theoretical calculations 
was considered excellent for flows involving pressure disturbances 
(including compression shocks or expansion zones), temperature— contact 
discontinuities, and reflections or interactions of these quantities. 
Slight discrepancies were observed in the correlation of the data and were 
believed to be introduced by failure to generate exactly the idealized 
type of flow assumed in the theory. The successful recording of the 
unsteady flow properties was largely attributed to the special equipment 
developed for this purpose. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., March 17, 19^9 
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APPENDIX 

•r 

INSTRUMENTATION 

Apparatus for Measuring and Recording Pressures 


Crystal, preamplifier, and filter .— The theoretical shape of the 
pressure history to he recorded is shown in figure 10. The steepness of 
the leading edge and the flatness of the constant— pressure portion of 
the record are dependent upon the response characteristics of the crystal 
and preamplifier combination. The crystal and preamplifier combination 
required to reproduce the pressure history shown is essentially one with 
a very low -frequency response as well as a high-frequency response. 

Tne pressure-measuring element used was a Massa M— 101 piezoelectric 
pressure pick-up . The first records taken showed that the high-frequency 
response of the crystal was sufficient to reproduce satisfactorily the 
steep leading edge but that poor low-frequency response resulted in a 
noticeable drop in the output during the time of constant pressure. The 
open— grid cathode— follower preamplifier circuit shown - in figure 15 was 
developed to present a maximum impedance to the crystal unit in order to 
obtain low-frequency response. The addition of the 100-mi cromi cr of arad 
condenser across the input further increased the low-frequency response 
without noticeably affecting the steepness of the wave front. 

Figure 1 6 shows the input— output voltage of the preamplifier alone. 
Negative input voltages are seen to have resulted in distortion which 
made it necessary to bias the grid positively when recording calibrating 
voltages if a sine-wave voltage was used. The grid- to- cathode resistance 
is shown to be effectively infinite until the grid is 14 volts positive 
so that leakage of the piezoelectric charge was due only to the internal 
resistance of the crystal itself. Records of crystal leakage were made 
by impressing a direct— current voltage momentarily across the crystal 
when in combination with the preamplifier and 100-micromicrofarad con- 
denser and by recording the exponential decay of the voltage output of 
the preamplifier; the procedure was then repeated without the crystal in 
the circuit. The time constant of the crystal was found to be 55 milli- 
seconds; that is, two— thirds of the charge on the crystal was dissipated 
in 55 milliseconds, about 1;=;- percent per millisecond. In order to maintain 
this time constant, the crystal must be dehydrated after each day's use 
by storing it in a flask containing a dessicant. The slope of the trace 
during the constant— pressure portion of the records shown in figure 10 
illustrates the very sma l l leakage occurring during this time. 

The appreciable amplitude of hash observed on the record shown in 
figure 17 (a) was generated by the crystal's resonating at its natural 
frequency of about 47 kilocycles after being excited by the steep wave 
front. Satisfactory attenuation of the resonant frequency was obtained 
by the use 6f filters as is shown by the record of figure 17(t>) which 
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was taken under conditions similar to those for figure 17(a) except that 
a filter was used. The filter was constructed in accordance with the 
diagram shown in figure l8 with parts that were accurate to within ±1 per- 
cent and functioned satisfactorily with the different crystals used 
during the investigation. With the filter in the circuit some attenuation 
of frequency components above 5000 cycles can be expected; however, the 
rise time at the leading edge would be the only part affected in these 
records . 

Because the preamplifier output for the range of pressures covered in 
these tests was of the order of 10 volts, amplification, in addition to 
that of the oscilloscope amplifier, was not needed. 

Becording apparatus .— At first, records were made by photographing 
a four— beam commercially constructed oscilloscope. The four beams made 

it possible to record simultaneously the pressure history, a — second 

timer, a 1000 cycle calibrating voltage, and a reference line. Because 
defocusing at high intensities masked many details in the records made 
with the four— beam oscilloscope, it was replaced with a single— beam 
oscilloscope of higher resolution that gave records of more detail and 
higher contrast. When this single— beam oscilloscope was used, the 
reference line and calibrating voltages were recorded immediately preceding 
a test on the same film record. Calibration indicated that the spot 
deflection on the oscilloscope varied linearly with input voltage to the 
preamplifier until saturation of the preamplifier. 

Moving film and no sweep of the spot, rather than single sweep with 
stationary film, was chosen as the method of photographing the record. 
Eecords were taken with an NACA drum camera using an f/ 2.0 lens and 

orthochromatic film 2i — inches wide. This camera consists of a drum 

lo 

l6j£ inches in circumference rotated at 3^00 rpm by means of a synchronous 

motor so that the film speed is 0.95 ±- 1 percent inch per millisecond. 

The shutter was set for l /50 second and was synchronized with the drum 
of the camera and solenoid trigger on the experl mental apparatus as 
described in the following section. 

Synchronizing circuit .— Figure 19 shows a block diagram of the 
synchronizing circuit. A small bar magnet cemented on the film drum 
generated a trigger signal as it passed a coil mounted on the drum-case 
lid. The trigger signal caused the shutter solenoid to be energized 
which in turn released the shutter. After a delay of a few milliseconds, 
the external solenoid for actuating the diaphragm-rupturing trigger on 
the experimental apparatus was also energized. Film was loaded on the 
drum so that it began and ended at the bar magnet, thereby establishing 
a definite time relation between the trigger signal and the instant the 
beginning of the record passed the lens. The shutter was set for 
l/50 second or one complete revolution of the film drum. 
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The pulse-delay circuit (fig- 20 ) which synchronized the schlieren 
apparatus with the diaphragm "burst was actuated by a crystal near the ' * 

diaphragm and was set to give the delay required for the shock to travel 
from the crystal to the glass section of the experimental apparatus. 

The delay was set by means of a 10 — turn helical potentiometer which was 
calibrated by use of a Potter interval timer. It was determined that 
the delay interval could be reset within ±10 microseconds consistently . 

The pulse— delay circuit was connected to the socket of existing spark 
equipment used for schlieren photography so that the output voltage 
pulse fired a thyratron which was in series with the primary winding of 
an ignition coil. The ignition coil in turn fired a series pilot gap 
which -set off the main spark gap. 


Apparatus for Recording Shock Speed 

The speed of the shock was determined by measuring the interval 
between the time the shock passed the window of the experimental apparatus 
and the time it reached the end plate by means of the following arrange- 
ment: A sine wave of 10,000 cycles per second was impressed on the 

vertical axis of an oscilloscope to provide a time scale, and two pulses 
representing these two positions of the shock were applied to the 

horizontal scale. One of these pulses was the variation in the output * - 

of a 931 photomultiplier tube caused by the bending by the shock wave of 
a beam of light which emerged from a slit 0 . 020 — inch wide and 1 — inch 
high, passed through the windows of the experimental apparatus at a » 

point 19 ^ inches from the end of the low-pressure chamber, and activated 

the photomultiplier tube. (The light from the slit was cut off from the 
photomultiplier by a knife-edge when no Bhock was present.) The other 
pulse was the output of a crystal caused by the impact of the shock on 

the end plate of the 72 ^- -inch low-pressure chamber to which the crystal 

was attached. The trace on the oscilloscope was photographed with the 
drum camera synchronized to the diaphragm burst as previously described. 

The sine-wave frequency . was accurate to within one part in 100,000 and 
the outputs of the photocircuit and crystal (without preamplifier and 
filter) were sharp enough for the departure of the spot from the sine 
pattern to be read to about 1 microsecond in 1200. Six records were 
taken at each point shown in figure 2 and the maximum scatter of the 
records at each point was less than ±1 microsecond in about 1200 . 

I 
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TABLE I 

THEORETICAL CALCULATIONS OF UNSTEAUT FLOW VALUES 
(Calculated for T Q = ^ = 67. V 3 f] 


U£ 

a o 

p 2 

p o 

CVJ 

CVJ 

1 

CVJ 

3 


*2 

p l 

p 2 “ p l 

T 0 ” T 22 
(°F ) 

T 2 -Ti 
(°F ) 

Yo 

h 

P 2 

h 

p 0 

P3 - P2 

t 3“ T 2 

(°F) 

a o 

aQ 

Pi 

p 0 

p 0 

a 0 

p o 

a 0.10 

0.868 

-0.880 

1.062 

1.149 

0.756 

0.112 

20.9 

21.3 

-0.981 

1.147 

0.994 

0.126 

21.8 

.1445 

.816 

-.828 

1.089 

1.223 

.66 7 

.148 

29.8 

30.7 

-.9 755 

1.216 

.989 

.1755 

32.0 

.15 

.808 

-.820 

1.094 

1.230 

.657 

.151 

31.2 

32 ;2 

-.974 

1.222 

.988 

.180 

33.2 

a .20 

.751 

-.760 

1.127 

1.316 

.571 

.180 

41.4 

43.3 

-.966 

1.302 

.979 

.227 

45.0 

.248 

.701 

-.703 

1.159 

1.402 

.500 

.201 

51.0 

54.2 

-.9612 

1.380 

.968 

.266 

57.1 

.25 

.69 8 

-.700 

l.l6l 

1.4o6 

.497 

.202 

51.5 

54.6 

-.961 

1.384 

.967 

.268 

57.3 

a .30 

.649 

-.640 

1.196 

1.502 

.432 

.217 

61.4 

66.2 

-.955 

1.469 

.953 

.304 

70.0 


Values taken from reference 8. 
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Uq = ° 


U 2 = 0 

(a) t = 0. 




Figure 1.— Pressure distribution at typical times along constant area 
passages which were separated hy diaphragm at t = 0. 


Pressure ratio across shock, pg /p^ 


Value of P 2 /p^ calculated from 
equation (5) and measurements of 
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Pressure ratio across diaphragm, Pq/p^ 


Figure 2 .— Theoretical and experimental pressure ratio across compression shock before reflections 
as function of pressure ratio across diaphragm before bursting. &q = a-^. 
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u = 0 { 
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p = 1.40 
a = 1070 

n — oon 



Ip 

iH 

iPittl 

U - ^oU 

— v 


Step 1.- Zone AB considered as zones Y and Z. 


Tube 

end 


p = 2.0 

a = 1126 
u = 0 


Front Y 


1048 


P = 1.70 
a = 1100 
u = 130 


Front Z 


880 


p = 1.40 
a = 1070 
u = 280 


Step 2.- Zones Y and Z replaced with fronts Y and Z. 


Tube 

end 



p = 1.445 Front Y 

Front Z 


a = 1074 reflected 



u = 0 < — 

y 


— t — * 



[1152 * 



Step 3.- After reflection of front Y. 


Tube 

end 


“1 p = 1.191 r 

Front Z a = 1044 Front Y 

<—\ — u = 150 (— ► > 

984 ^ 1272 


Step 4.- After interaction of reflected front Y and front Z. 


Tube 

end 


p = 0.968 
a = 1014 
u = 0 


Front Z 

reflected 

— h-> 

1104 


Front Y 

> 



Step 5.- After reflection of front Z. 


Figure 3.— Illustration of treatment of expansion zone "by considering 
the expansion zone as a number of discontinuous expansion fronts. 
The values shown were obtained by assuming the fluid to be air 
(7 = 1 . 4 ) . 
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Denotes these quantities constant 
throughout an enclosed area 



Figure h.— Typical diagram of x against t showing position histories 
of unsteady flow quantities following diaphragm hurst in constant 
area passages, ag = a^. 


Pressure 



Figure 5 . Pressure distribution along axis of propagation of expansion zone. 
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Figure 6.— Pressure variation with time at reflection of expansion zone. 


NACA TN I9O3 



1 

] X End plate 


(a) Compression and expansion chambers. 


Crystal output 
terminals 



Solenoid support 


Rubber seal 


Crystal face 


Side wall of 8 -inch channel section 
Diameter = g^-in. 

(c) Trigger and solenoid mount. 


PHP 


(b) Crystal unit mount. 

Figure J Sketch of apparatus used to study flow. 
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Figure 8.— Apparatus used to study flow 
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(a) 24^ — inch low-pressure chamber. 

Figure 9.- Distance-time diagram for various low-pressure chambers 
£0 = 2.0; a Q = a-L = 1126 feet per second. 


Time, milliseconds 
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(t>) 48 — inch low-pressure chamber. 
Figure 9.— Continued. 
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(c) 



— inch low-pressure chamber . 


Figure 9.— Concluded. 
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(a) Crystal located in end plate of 24 i -inch low-pressure chamber. 



Figure 10.— Theoretical and experimental pressure— time histories for 
various low-pressure chambers and crystal positions plotted on same 

pressure— time scales. — = 2.0; a Q = a 1 = 1126 feet per second. 
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Figure 10.— Concluded 
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(a) 


- ij = 2.0: 48— inch low— pressure chamber at end plate; end plate 

p i - 

- 2 -— inch thick plus 6 -pound weight. 





(d) ^ = 2 . 5 ; 24^— inch low-pressure chamber at end plate. 

Figure 11.— Typical pressure records obtained with standard test configu- 

ration except where noted. 
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(e) 


— = 1.3; 24“ — inch low-pressure chamber at end plate. 
P 1 2 



side 24 inches from diaphragm in low-pressure chamber. 



side 24 inches from diaphragm in low-pressure chamber. 
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2.0; 72jj — Inch low-pressure chamber; crystal located in the 
side 20 inches from diaphragm in high-presaure chamber. 



2.0; 72^— inch low-pressure chamber; crystal located at end 
plate of high-pressure chamber. 



1^3; 72I— inch low-pressure chamber; crystal located at end 


plate of high-pressure chamber. 
Figure 11 .— Concluded. 
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INFORMATION REGARDING TEST CONDITIONS ON CURVES 



Crystal output, volts 
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O .2 .4 .<5 .8 /.O /.2 /.4 /.<5 

Ap/p a (calculated from Po/p^) 


Figure 12.— Voltage output of crystal unit as a function of pressure 

difference experienced at crystal as a result of passing disturbance. 
Ap calculated from theory. (Symbols designating test conditions 
are identified in table II.) 



Zero voltage line 



Figure 13.— Method used to measure voltage on pressure— time records. 
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(a) — = 1.6; 2^=— inch low-pressure chamber; 

first window from diaphragm; time, 

0.5 millisecond. 



(h) ^-2 = 2.0; 24— — inch low-pressure chamber; 

Pi 2 

first window from diaphragm; time, 

0.5 millisecond. 



(c) — = 2.0; 2ti — inch low-pressure chamber; 

Pi 2 

second window from diaphragm; time, 

0.97 millisecond. 


Figure 14.— Schlieren photographs taken at various times following "burst 
of diaphragm. Positive x— direction is to the right . 


^ L-5981J.6 





NACA TN I9O3 


55 



(d) — = 2.0; 2 h~ — inch low-pressure chamber; 

Pi 2 

second window from diaphragm; time, 

1.08 milliseconds. 



(e) — = 2.0; 2 b~- inch low-pressure chamber; 

P 1 2 

second window from diaphragm; time, 

I.98 milliseconds; shock reflected. 



(f) — = 2.0; 2 ^ — inch low-pressure chamber; 

Pi 2 

first window from diaphragm; time, 

2.0 milliseconds. 

Figure 1^.— Continued. 



l-5981;7 




NACA TN I9O3 


57 



(g) 


— = 2.0; 2 4 ^— inch low— pressure chamber; 

Pi 2 


first window from diaphragm; time, 
2.7 milliseconds; shock reflected. 


00 



— = 2.0; 72^— inch low-pressure chamber; 

Pi 2 

second window from diaphragm; time, 

1.22 milliseconds. 



(i 1 ! = 2.0; 72— —inch low-pressure chamber; 

Pi 2 

second window from diaphragm; time, 

5.3 milliseconds; shock reflected. 


Figure ll.— Concluded. 
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Figure 15.— Preamplifier circuit. 
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Figure 16.- Input-output voltage characteristics of the preamplifier. Plate supply voltage 

90 volts; cathode resistance, 2200 ohms. 
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(a) Record taken without filter in circuit. 



(b) Record taken with filter in circuit. 


Figure 17.— Experimental records showing attenuation of natural frequency 

of crystal by 4 7 - kilocycle filter. 
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Figure 18.- Filter for removing the crystal frequency of ^7 kilocycles 
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drum lid 
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Closing switch applies positive peak to thyratron 
grid causing it to fire and operate relay (l). 

Rela;y (l) connects the shutter-solenoid supply to the 
shutter solenoid to take the picture. Relay (l) also 
disconnects the main supply from the thyratron to prevent 
interference with preamplifier. Thyratron current is 
then obtained from the 8-microfarad condenser connected 
to its plate. 

Relay (2) operates as soon as the shutter 
solenoid supply voltage has dropped a sufficient 
amount below the main supply voltage and connects 
the external solenoid to its supply. Relay (2) 
also extinguishes the thyratron by opening its 
cathode line. 



High voltage 
to preamplifier 

i 


External 

solenoid 


External- 


Shutter- 

solenoid 


solenoid 

high- 


high- 

voltage 



voltage 

supply 



supply 


Resistor to prevent 
rapid discharge 
between solenoid 
supplies 


Figure 19.— High-speed camera synchronizing unit. 
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Figure 20.— Circuit for 180 to 1400 microsecond pulse delay. 

(K, 1000 ohms; w, watts; meg, megohms; }if, microfarads; v, volts.) 


NACA TN I 9 O 3 


4 


% 




